Measurement of ep — > e p7r + 7r and baryon resonance analysis 
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The cross section for the reaction ep — > e'piY + iv~ was measured in the resonance region for 
1.4<W<2.1 GeV and 0.5< Q 2 <1.5 GeV 2 /c 2 using the CLAS detector at Jefferson Laboratory. 
The data shows resonant structures not visible in previous experiments. The comparison of our 
data to a phenomenological prediction using available information on TV* and A states shows an 
evident discrepancy. A better description of the data is obtained either by a sizeable change of 
the properties of the Pi3(1720) resonance or by introducing a new baryon state, not reported in 
published analyses. 

PACS numbers: 13.60.Le, 13.40.Gp, 14.20.Gk 



Electromagnetic excitation of nucleon resonances is 
sensitive to the spin and spatial structure of the transi- 
tion, which in turn is connected to fundamental proper- 
ties of baryon structure, like spin-flavor symmetries, con- 
finement, and effective degrees of freedom. In the mass 
region above 1.6 GeV, many overlapping baryon states 
are present, and some of them are not well known; mea- 
surement of the transition form factors of these states 
is important for our understanding of the internal dy- 
namics of baryons. Many of these high-mass excited 
states tend to decouple from the single-meson channels 
and to decay predominantly into multi-pion channels, 
such as A7T or Np, leading to Ntttt final states [Q. 
Moreover, quark models with approximate (or "broken") 
SU(6)(g>0(3) symmetry fi || predict more states than 
have been found experimentally; QCD mixing effects 
could decouple these unobserved states from the pion- 
nucleon channel while strongly coupling them to two- 
pion channels ^, ]|, ||] . These states would therefore not 
be observable in reactions with irN in the initial or final 
state. Other models, with different symmetry properties 
and a reduced number of degrees of freedom, as e.g. in 
ref. ||, predict fewer states. Experimental searches for at 
least some of the "missing" states predicted by the sym- 
metric quark models, which are not predicted by models 
using alternative symmetries, are crucial in discriminat- 
ing between these models. Electromagnetic amplitudes 
for some missing states are predicted to be sizeable |2| as 
well. Therefore, exclusive double-pion electroproduction 
is a fundamental tool in measuring poorly known states 
and possibly observing new ones. 

In this paper we report a measurement of the ep — > 
e'p7r + 7r~ reaction studied with the CEBAF Large Ac- 
ceptance Spectrometer (CLAS) at Jefferson Lab. Beam 
currents of a few nA were delivered to Hall B on a liquid- 
hydrogen target, corresponding to luminosities up to 
4 x 10 33 cm _2 s _1 . Data were taken in 1999 for about two 



months at beam energies of 2.6 and 4.2 GeV. Important 
features of the CLAS § are its large kinematic coverage 
for multi-charged-particle final states and its good mo- 
mentum resolution (Ap/p ~1%). Using an inclusive elec- 
tron trigger based on a coincidence between the forward 
electromagnetic shower calorimeter and the gas Cercnkov 
detector, many exclusive hadronic final states were mea- 
sured simultaneously. Scattered electrons were identi- 
fied through cuts on the calorimeter energy loss and the 
Cerenkov photo-electron distribution. Different channels 
were separated through particle identification using time- 
of- flight information and other kinematic cuts. We used 
the missing-mass technique, requiring detection in CLAS 
of at least epir + . The good resolution allowed selection 
of the exclusive final state, epTr + ir~ . After applying all 
cuts, our data sample included about 2 x 10 5 two-pion 
events. 

The range of invariant hadronic center-of-mass (CM) 
energy W (in 25 MeV bins) was 1.4-1.9 GeV for the first 
two bins in the invariant momentum transfer Q 2 , 0.5- 
0.8 (GeV/c) 2 and 0.8-1.1 (GeV/c) 2 , and 1.4-2.1 GeV 
for the highest Q 2 bin, 1.1-1.5 (GeV/c) 2 . Data were 
corrected for acceptance, reconstruction efficiency, radia- 
tive effects, and empty target counts. They were further 
binned in the following set of hadronic CM variables: in- 
variant mass of the pir + pair (10 bins), invariant mass 
of the 7r + 7T~ pair (10 bins), ir~ polar angle 6 (10 bins), 
azimuthal angle <f> (5 bins), and rotation freedom tf> of 
the pir + pair with respect to the hadronic plane (5 bins). 
The full differential cross section is of the form: 



da 

dWdQ 2 dAI p7T +dM 7T + 7r - dcos 0^- d<f> v - dip p7r + 



dM pn + dM 7r+7r - d cos 6^-d<j) v - dyb pyr+ 
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r„ = 



a 1 W(W 2 - Ml) 
4^E 2 M 2 (l-e)Q 2 



(2) 



where r„ is the virtual photon flux, is the virtual pho- 
ton cross section, a is the fine structure constant, E is 
the electron beam energy, M p is the proton mass, and e is 
the virtual photon transverse polarization || . Systematic 

30 
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FIG. 1: Left: Total cross section for "f v p ~ * pix + -K~ as a 
function of W. Data from CLAS are shown at Q 2 =0.5- 
0.8 (GeV/c) 2 (full points), Q 2 =0.8-l.l (GeV/c) 2 (open 
squares), and Q =1.1-1.5 (GeV/c) 2 (open triangles). Error 
bars are statistical only, while the bottom band shows the 
sytematic error for the lowest Q 2 bin. The curves represent 

^— from 



our step (A) reference calculations. Right: 

CLAS at Q 2 =0.8-l.l (GeV/c) 2 and W=l. 7-1. 725" GeV (sta- 
tistical error bars only). The curves represent our step (A) 
reference calculations, extrapolated to the edge points. The 
dashed line includes all resonances, the dot-dashed line in- 
cludes only the non-resonant part, and the solid line is the 
full calculation. 

uncertainties were estimated as a function of W and Q 2 . 
The main sources were acceptance modeling, finite inte- 
gration steps, and modeling of the radiative corrections, 
each one being at the 3 to 10% level. Each of the various 
cuts applied (fiducial, missing mass, etc.) contributed 2 
to 5%. In Fig. [I] (left) we report the total virtual photon 
cross section as a function of W for all Q 2 intervals ana- 
lyzed. The CLAS data points clearly exhibit structures, 
not visible in previous data H due to limited statistical 
accuracy. 

Since existing theoretical models jl(| are limited to 
W <1.6 GeV, we have employed a phenomenological 
calculation jllj] for a first interpretation of the data. 
This model describes the reaction j v p — > pir + ir~ in the 
kinematic range of interest as a sum of amplitudes for 



IvP 



Air — > pir + ir and "f v p 



P7T + 7T 



while 



all other possible mechanisms arc parameterized as phase 
space. A detailed treatment was developed for the non- 
resonant contributions to Air, while for pp production 
they were described through a simple diffractive ansatz. 
For the resonant part, a total of 12 states, classified as 
4* fjj], with sizeable Aw and/or pp decays, were included 



based on a Breit-Wigner ansatz. A few model parameters 
in non-resonant production were fitted to CLAS data at 
high W , where the non-resonant part creates a forward 
peaking in the angular distributions, and kept fixed in 
the subsequent analysis. The phase between resonant 
and non-resonant Air mechanisms was fitted to the CLAS 
data as well. To simplify the fits, we reduced eqn. (1) to 
three single-differential cross sections, dA d<T - , dM da + — , 

" T , by integrating over the other hadronic vari- 



and 



d cos 9 

ables. These three 1-D distributions were then fitted si- 



multaneously. Here 



dAI 



and 



da 



I COS P _ 



are both con- 



nected with the dominant A ++ ir production reaction, 
while dM d<T + — is connected with pp° production. For 

each W and Q 2 bin, a total of 26 data points from the 
three single-differential cross sections were used in our 
fits. The two edge points in both the pir + and ir + ir~ 
mass distributions were excluded as the model did not 
take into account the kinematic smearing in the M„+ n - 
versus M„„+ plot caused by the W bin width. 

The data analysis was performed in the following steps: 
(A) We produced reference curves using the available in- 
formation on the N* and A resonances in 1.2-2 GeV mass 
range. Discrepancies between the CLAS data and our 
calculation were observed, which led to the subsequent 
steps B and C. (B) Data around W=l.7 GeV were fit- 
ted using the known resonances in the PDG but allowing 
the resonance parameters to vary in a number of ways. 
The best fit, corresponding to a prominent Pi 3 partial 
wave, could be attributed to the PDG Pi3(1720) reso- 
nance, but with parameters significantly modified from 
the PDG values. (C) As an alternative description, we 
introduced a new baryon state around 1.7 GeV. In what 
follows we describe each of the steps above in more detail. 

evo- 



Step (A) - To produce our reference curves, the Q 



lution of the A1/2 and A 3 / 2 electromagnetic couplings for 
the states was taken either from parameterizations of ex- 
isting data ]12|], or from Single Quark Transition Model 
(SQTM) fits Jl2| where no data was available. For the 
Pti(1440) (Roper), given the scarce available data, the 
amplitude A 1 / 2 was taken from a Non-Relativistic Quark 
Model (NRQM) Q. Partial LS decay widths were taken 
from a previous analysis of hadronic data J LJ] and renor- 
malized to the total widths from Rcf. p|. Results for 
step (A) are reported in Fig. |l| (left). The strength of 
the Di3(1520) resonance at 1.5 GeV and the underlying 
continuum are well reproduced, except for the region on 
the low-mass side of the peak. However, a strong discrep- 
ancy is evident at W around 1.7 GeV. Moreover, at this 
energy the reference curve exhibits a strong peak in the 
ir + ir~ invariant mass (Fig. [l], right), connected to size- 
able p meson production. This contribution was traced 
back to the 70-85% branching ratio of the Pi3(1720) into 
this channel |IJ [ij, HJ. 

Step (B) - Starting from the above mentioned reference 
values, the parameters of various states were varied in or- 
der to fit the CLAS data. In this discussion, we restrict 
ourselves to the discrepancy around 1.7 GeV and the few 
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TABLE I: PDG Pi 3 (1720) parameters from fit (B) and new 
state parameters from fit (C). Errors are statistical. 





M (MeV) 


T (MeV) 


^ (%) 


£f£ (%) 


PDG P13 (B) 


f725±20 


ii4±i9 


63±12 


19±9 


PDG f] 


f 650- f 750 


100-200 


N/A 


70-85 


new P13 (C) 


f720±20 


88±17 


41±13 


17±10 



resonant states relevant in this energy region. All fit x 2 jv 
values were calculated from the 8 W bins between 1.64 
and 1.81 GeV and from the 3 Q 2 bins (624 data points). 
The number of free parameters ranged from 11 to 32, de- 
pending on the fit, corresponding to ^=613 to 592 degrees 
of freedom. Assuming the resonance properties given by 
the PDG, the bump at about W—1.7 GeV cannot be 
due to the L»i 5 (1675), Pi 5 (1680), or L> 33 (1700) states; 
the first because its well known position cannot match 
the peak; the second because of its well known position 
and photocouplings Jl6| ; the third due to its large width 
(^300 MeV). The remaining possibilities from the PDG 
are the L>i 3 (1700), the P 13 (1720), and the P u (1710) (the 
latter not present in step (A)), for which no data on Ai/ 2 
or at Q 2 > are available [|l6|. If no configuration 
mixing occurs, the Pi 3 (1700) cannot be excited in the 
SQTM, while the SQTM prediction for the Pi 3 (1720) re- 
lies on ad hoc assumptions. According to the literature 
jl], |TJ, jl5| , hadronic couplings of the £>i 3 (1700) and the 
total width of the Pn(1710) are poorly known, while the 
Pi 3 (1720) hadronic parameters should be better estab- 
lished. Several other partial waves were investigated in 
step (C). 




ried out: the Pn(1440) strength was fitted to our low W 
data; the £>is(1675) and the £>i 3 (1700) photocouplings 
(which vanish in the SQTM) were replaced by NRQM 
values from Ref. |ll|; an empirically established ^1/2, 3/2 
SQTM fitting uncertainty or NRQM uncertainty of 10 
or 20% (a) was applied to all N* states; the hadronic 
parameters were allowed to vary for the Di 3 (1700) ac- 
cording to Ref. |l4|] ; and finally, the curves providing the 
best x 2 /^ were selected as the starting points. Then, first 
we performed three fits, (Bl), (B2), and (B3), where the 
photocouplings of only one resonance at a time were var- 
ied. In (Bl), we varied Ax/2, A 3 / 2 , hadronic couplings, 
and position of the Di 3 (1700) in a wide range. In (B2), 
the same was done for the Pi 3 (1720), and in (B3) for the 
Pn(1710). In both fits (B2) and (B3), we also varied the 
hadronic parameters and the position of the Di 3 (1700) 
over a range consistent with their large uncertainties from 
Ref. jL4j. Fits (Bl) and (B3) gave a poor description 
of the data, with x 2 /z/=5.2 and 4.3, respectively. The 
best fit { X 2 /v =3.4) was obtained in (B2) (Fig.|). How- 
ever, the resulting values for the branching fractions of 
the P 13 (1720) were significantly different from previous 
analyses reported in the literature and well outside the re- 
ported errors 

§, |§ H. Starting from (B3), we then per- 
formed a final fit, (B4), for which the Pi 3 (1720) hadronic 
couplings were fixed from the literature, and varying the 
photocouplings of all three candidate states, Di 3 (1700), 
Pi 3 (1720), and Pn(1710), by 100% (a). No better solu- 
tion was found, the \ 2 l v being 4.3 (Fig. ||). In Fig. || we 
report the final comparison of fits (B2) and (B4) with the 
total cross section data. Table | shows our results (first 




to bottom) at W=\. 7-1. 725 GeV and for the three mentioned 
Q 2 intervals (left to right). The error bars include statisti- 
cal errors only. Curves (see text) correspond to the fits (B2) 
(red), and (B4) (blue), and are extrapolated to the mass dis- 
tributions edge points. 

To improve our reference curves before fitting the 
bump at around 1.7 GeV, the following steps were car- 



1.6 1.8 

W (GeV) 

FIG. 3: Left: Total cross section for j v p — > pn + ir~ as a 
function of W from CLAS at the 3 mentioned Q 2 intervals 
(see fig.l). The error bars are statistical only. The curves 
CN?(deg) ( see text) correspond to the fits (B2) (red), and (B4) (blue). 

Right: subdivision of the fitted cross section (B2) for Q 2 =0.5- 
from CLAS (from top q g (QeV/c) 2 into resonant A++TT' (black solid), continuum 



A ++ 7r~ (black dashed), resonant p°p (magenta solid), and 
continuum p°p (magenta dashed). Notice the different verti- 
cal scales. 

row) with statistical uncertainties, in comparison with 
the PDG values (second row). Our fits were not pro- 
viding an unambigous separation of Ax/2, A 3 / 2 i an d the 
longitudinal Si/ 2 , so we report as result the total pho- 
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TABLE II: PDG Pi 3 (1720) total photocoupling from fit (B2) 
and new state total photocoupling from fit (C). Errors are 
statistical. 



step 


Q 2 

(GeV/c) 2 


\Z^l/2 + ^3/2 + ^1/2 

(10- 3 /VGeV) 


B2 


0.65 


83±5 


B2 


0.95 


63±8 


B2 


1.30 


45±27 


C 


0.65 


76±9 


C 


0.95 


54±7 


c 


1.30 


41±18 



tocoupling strength, JA^^ + A^ 2 + S^ T Such value 

for the Pi3(1720) fit is reported in the first three rows 
of Table |J. The errors reflect the statistical uncertain- 
ties in the data and the correlations among the different 
resonances. 

As discussed above, fitting the data around 1.7 GeV 
with established baryon states leads either to a poor fit or 
to a drastic change in resonance parameters with respect 
to published results. In the framework of our analysis, 
there is no way to assess the reliability of the previously 
determined hadronic parameters of the PDG Pi3(1720). 
The resonant content of the reaction ttN — * irirN, which 
is used to obtain the hadronic parameters, may be dif- 
ferent from that of reactions initiated by an electromag- 
netic probe. In particular, the Pi3(1720) state seen in 
ttN — > ttttN may not be excited in electroproduction, 
while some other state that decouples from ttN may be 
excited electromagnetically. This possibility is studied in 
the next step. 

Step (C) - We investigated whether our data could be 
fitted by including another baryon state, while keeping 
the hadronic parameters of the Pi3(1720) as in Refs. 
Jj], The quantum numbers Sji,Pn,Pi3,Dj^,Dj5, 

Fi5,Fn were tested on an equal footing, where 7/2 is the 
isospin, undetermined in our measurement. We then si- 
multaneously varied the photocouplings and the hadronic 
parameters of the new state and the -Di3(1700). The to- 
tal decay width of the new state was varied in the range of 
40-600 MeV, while its position was varied from 1.68-1.76 
GeV. The best fit (x 2 /i/=3.3) was obtained with a P/3 
state, while keeping the Pi 3 (1720) hadronic parameters 
at published values. Other partial waves gave a^ 2 /" > 
4.2. Curves obtained from the best fit were nearly iden- 
tical with the red solid lines in Figs. ^ and |3|. In order to 
avoid the unobserved p production peak (Fig. ^, right), 
the photocouplings of the PDG Pi3(1720) had to be re- 
duced by about a factor of two with respect to the SQTM 
prediction, making its contribution very small. Instead, 
in this fit the main contribution to the bump came from 
the new state. Resonance parameters and total photo- 



coupling value obtained for the hypothetical new state 
are reported in Table | (last row) and [H] (last 3 rows), 
respectively. 

A second P13 state is indeed predicted in the quark 
model of Ref. Q], however with a mass of 1870 GeV. The 
presence of a new three-quark state with the same quan- 
tum numbers as the conventional Pi3(1720) in the same 
mass range would likely lead to strong mixing. How- 
ever, as mentioned above, a different isospin and/or par- 
tial wave cannot be excluded. Finally, the new state 
may have a different internal structure, such as a hybrid 
baryon with excited glue components. Such a P13 hybrid 
state is predicted in the flux tube model JlTj . Yet another 
possibility is that some resonance parameters established 
in previous analyses may have much larger uncertainties 
than reported in the literature. In this case, outlined 
in our step (B), our analysis would establish new, more 
precise parameters for a known state, and invalidate pre- 
vious results. 

In conclusion, in this letter we presented data on the 
ep — > e'p7r + 7r~ reaction in a wide kinematic range, with 
higher quality than any previous double pion produc- 
tion experiment. Our phenomenological calculations us- 
ing existing PDG parameters provided a poor agreement 
with the new data at W ~ 1700 MeV. We explored 
two alternative interpretations of the data. If we dis- 
miss previously established hadronic parameters for the 
Pi3(1720) we can fit the data with a state having the 
same spin/parity/isospin but strongly different hadronic 
couplings from the PDG state. If, alternatively, we in- 
troduce a new state in addition to the PDG state with 
about the same mass, spin |, and positive parity, a good 
fit is obtained for a state having a rather narrow width, a 
strong An coupling, and a small pN coupling, while keep- 
ing the PDG Pi3(1720) hadronic parameters at published 
values. In either case we determined the total photocou- 
pling at Q 2 > 0. A simultaneous analysis of single and 
double-pion processes provides more constraints and may 
help discriminate better between alternative interpreta- 
tions of the observed resonance structure in the CLAS 
data. Such an effort is currently underway. 
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